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We explore the nature of the faint blue objects in the Hubble Deep Field 
South. We have derived proper motions for the point sources in the Hubble Deep 
Field South using a 3-year baseline. Combining our proper motion measurements 
C\J ■ with spectral energy distribution fitting enabled us to identify 4 quasars and 42 

stars, including 3 white dwarf candidates. Two of these white dwarf candidates, 
HDFS 1444 and 895, are found to display significant proper motion, 21.1 ± 7.9 
Q ' mas yr~^ and 34.9 ± 8.0 mas yr^^, and are consistent with being thick disk or 

^ ■ halo white dwarfs located at ~2 kpc. The other faint blue objects analyzed by 

■ Mendez & Minniti do not show any significant proper motion and are inconsistent 

with being halo white dwarfs; they do not contribute to the Galactic dark matter. 
O ■ The observed population of stars and white dwarfs is consistent with standard 

• Galactic models. 

^ ■ Subject headings: dark matter — Galaxy: halo — stars: evolution — white 

^ ■ dwarfs 

> 

. - - . 1. Introduction 

Faint blue objects discovered in deep Hubble Space Telescope images have been the 
subject of discussion in recent years. The extreme depth of the Hubble Deep Field (HDF) 
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North (Williams et al. 1996) and South (Casertano et al. 2000) and the Hubble Ultra Deep 
Field (HUDF; Beckwith et al. 2004) enables us to study faint stellar objects in the regions 
of the color-magnitude diagram that are devoid of standard Galactic stars. Supported by 
the observed microlensing events toward the Large Magellanic Cloud (Alcock et al. 2000), 
several investigators proposed that the faint blue objects observed in the HDF images can 
explain part of the dark matter in the Galaxy and a significant population of the halo of the 
Galaxy may be in the form of low-luminosity white dwarfs. Claims by Mendez & Minniti 
(2000; hereafter M&M) that the faint blue sources in the HDF North and South are Galactic 
stars seemed to be consistent with earher findings of Ibata et al. (1999) who found 5 faint 
halo white dwarf candidates with detectable proper motions in the HDF North. However, 
further analysis by Richer et al. (2001) and Kilic et al. (2004) showed that the faint blue 
objects in the HDF North do not show any significant proper motion. 

A detailed analysis of the point sources in the HDF North (Kihc et al. 2004) and 
the HUDF (Pirzkal et al. 2005) showed that blue extra-galactic sources may be confused 
with white dwarfs. The (lOcr) limiting AB magnitudes of the / band images for these 
two fields are 27.6 and 29.0, respectively. Using a 7 year baseline, Kilic et al. (2004) 
obtained proper motion measurements for the point sources in the HDF North including the 

5 faint blue objects; they identified two possible disk white dwarfs, one of which now also 
appears spectroscopically to be a white dwarf (D. Stern, private communication). Using low 
resolution spectroscopy and proper motion measurements, Pirzkal et al. (2005) identified 
20 late type stars, 2 quasars, and 4 possible white dwarf candidates in the HUDF. Kilic, 
von Hippel & Winget (2005) showed that only two of these candidates (HUDF 4839 and 
9020) are firm white dwarf candidates with HUDF 4839 possibly being a thick disk object, 
and HUDF 9020 being either a disk or a halo object. None of these white dwarf candidates 
show significant proper motion (Pirzkal et al. 2005). Non-detection of high velocity white 
dwarfs in the HUDF is consistent with non-dctcction of halo white dwarfs in the HDF North. 
Lack of detection of high velocity white dwarfs in these two fields implies that white dwarfs 
account for less than 10% of the Galactic dark matter (Pirzkal et al. 2005; Kihc, von Hippel 

6 Winget 2005). 

Ibata et al.'s (2000) and Oppenheimer et al.'s (2001) discoveries of apparent halo white 
dwarfs from kinematic surveys were enough to explain 2% of the dark matter in the solar 
neighbourhood. On the other hand, further analysis by several investigators showed that 
most of these white dwarfs are associated with the thick disk population of the Galaxy (Reid 
et al. 2001; Reyle et al. 2001; Bergeron et al. 2005). Nevertheless, old halo white dwarfs are 
observed in the Globular clusters M4 (Hansen et al. 2004) and NGC 6397 (Mendez 2002) 
and in the field toward M4 (Kalirai et al. 2004). Reid (2004; see for a complete review on 
high velocity white dwarfs) found that these observations do not require additions to the 
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standard Galactic populations. 

The HDF South data provides another opportunity to test whether the faint blue objects 
in deep Hubble images can be old halo white dwarfs and if they can explain part of the 
Galactic dark matter. M&M found 22 Galactic stars and 10 faint blue objects in the HDF 
South. If these 10 faint blue objects are halo white dwarfs, then they would explain 30-50% 
of the dark matter in the solar neighborhood. Wc extend Kilic et al.'s (2004) work to the 
HDF South by using the original HDF South data and images of the same field taken 3 years 
later for the GO-9267 proposal (WFPC2 Supernova Search, PI: S. Beckwith) to measure the 
proper motions of the point sources in the HDF South. Section 2 describes our first epoch 
data and the classification of the point sources, while data reduction procedures for our 
second epoch images are discussed in §3. §4 describes the proper motion measurements 
for the point sources. We present our spectral energy distribution fitting results in §5, and 
various implications of these results are then discussed in §6. 

2. Description of the Data 

The HDF South was imaged using the Wide Field Planetary Camera 2 (WFPC2) and 
the F300W (U300), F450W (B450), F606W (Veoe), and F814W (Igu) filters in 1998, and 
imaged again in 2001 using the WFPC2 and the F814W filter in order to search for high 
redshift supernova. The original images taken in 1998 reach down to AB = 26.8, 27.7, 28.3, 
and 27.7 in the four bands, respectively (lOcr, point source; Casertano et al. 2000). The 
(version 2) source catalogs for the first epoch are produced by the Space Telescope Science 
Institute (STScI) from the combined and drizzled images. Source detection was carried 
out on the inverse- variance- weighted sum of the F606W and F814W drizzled images. The 
combined F606W + F814W (hereafter V + I) image is significantly deeper than any of the 
combined images (Casertano et al. 2000). 

2.1. Identification of Point Sources 

M&M used the stellarity index from Source Extractor (SExtractor; Bertin & Arnout 
1996) to identify point sources and claimed that the star/galaxy separation is reliable for 
objects brighter than V + I =29. They visually inspected the objects with stellarity index 
<0.85 and found that they are extended. M&M identified all objects with stellarity index 
>0.90 as stars, and found 98 point sources in 4.062 arcmin^ of the HDF South. They analyzed 
objects with I < 27 (15(7 detections) and found 22 late type stars and 10 faint blue objects 
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(see M&M for a complete discussion). 

Kilic, von Hippel, and Winget (2005) showed that the SExtractor stellarity index fails 
at faint magnitudes. Instead, they suggested the use of the half light radius (R50) and the 
full width at half maximum (FWHM) measurements along with the empirical point spread 
function (PSF) distributions to identify the unresolved objects in deep HST images. Figure 

1 shows the distribiition of R50 with V + I magnitude for all HDFS objects (top panel) 
and for the objects with the stellarity index >0.85. Objects brighter than V + I = 20 are 
saturated and therefore have higher R50 and lower stellarity values. It is clear from this 
figure that the point sources arc well separated from the extended objects in this diagram 
down to V + I =28.5. A solid line (R50=3.0) marks the separation between the resolved 
and the unresolved objects. We note that we trimmed the source catalog to avoid spurious 
detections near the detector boundaries, and we also exclude the Planetary Camera (PC) 
observations because the PC has a different plate scale and the PSFs are different from the 
PSFs for the Wide Field Cameras. 

We use the IRAF task PRADPROF and the object centroids from Source Extractor 
to plot the radial profile of each object. Our empirical PSF distributions for objects with 
R50 < 3.0, the stellarity index > 0.85, and 21 < V + 1 < 28.5 are shown in Figure 2. As 
in M&M, we limit our analysis to the objects brighter than I = 27. The object number, 
R50, FWHM, and stellarity index of each object are shown in the top right corner of each 
panel. We also plot the PSF for a bright unsaturated, unresolved source (HDFS 10617) in 
each panel for a direct comparison. In addition, we examined the PSF distributions for the 
objects with R50 < 3.0 and the stellarity index < 0.85, and (as in M&M) found all of them 
to be resolved. 

A comparison of the PSF for each object with our template PSF (for HDFS 10617) shows 
that HDFS 126 (FWHM=6.03), 2003 (FWHM=5.38), and 497 (FWHM=5.28) arc clearly 
resolved. The typical FWHM for the unresolved objects is < 5. In addition, HDFS 1812, 
1827, 2129, and 567 have shallower PSF distributions than the template PSF, therefore 
are resolved, as well. Six of these seven objects have stellarity indices larger than 0.9, 
therefore would be classified as point sources based on the stellarity index. Our empirical 
PSF distribution analysis shows that the stellarity index from the SExtractor should not be 
used by itself to identify point sources in deep HST images. Instead, as Kilic, von Hippel, & 
Winget (2005) suggested, a combination of the R50, FWHM, stellarity index, and empirical 
PSF distributions is required to identify unresolved objects. This analysis leaves us with 
39 unresolved objects with 21 < V + I < 28.5, and / < 27, plus 7 brighter point sources 
{V + I< 21) in the HDF South. 

We note that we may have a few barely resolved objects in our sample, e.g. HDFS 
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191. Some extra galactic sources, e.g. Luminous Compact Blue Galaxies (LCBGs), can 
be quite small (half light radius=l-3 kpc) and blue {B — V < 0.6; Werk et al. 2004). 
Nevertheless, HST is able to resolve almost every galaxy if enough signal to noise is obtained. 
For example, Hoyos et al. (2004) could resolve four LCBGs at intermediate redshifts (z=0.1- 
0.44, R50 = 0.7-2.7 kpc) using WFPC2 {R50 = 0.23 - 0.47 "). In addition, using low 
resolution spectroscopy from the GRAPES survey, Pirzkal et al. (2005) showed that none 
of the unresolved objects in the HUDF (7 < 27) are LCBGs. Therefore, the number of 
compact galaxies that are classified as unresolved should be small. Even if our sample is 
contaminated by a few faint resolved objects, misclassification of these objects would not 
change the conclusions derived from our analysis (see §7). 

Figure 3 shows the distribution of the R50 and the FWHM for the HDF South objects 
in different magnitude ranges. Objects with stellarity index > 0.85 are shown as filled 
triangles, and the unresolved objects identified from their PSF distributions are shown as 
filled circles. All of the unresolved objects {V + I > 21) have stellarity index > 0.87, 
R50 < 3, and FWHM < 5. It is clear from this figure that the unresolved objects can 
be identified easily down to V" + / = 27. Even though, the classification becomes harder 
foT V + I > 27, five of the six objects that are fainter than V + I = 27 and classified as 
unresolved (HDFS 261, 441, 1020, 1306, and 2178) all have PSFs, R50, FWHM, and the 
stellarity values consistent with unresolved objects. Therefore, their identification as point 
sources is secure. The identification of HDFS 2488 (i?50 =2.96, I =26.97) as unresolved is 
questionable, but still included in our analysis for completeness purposes. Astrometric and 
photometric data for the 46 objects that we identified as unresolved are given in Table 1. 
We adopted the cahbrated photometry of Casertano et al. (2000). 

3. The Second Epoch Data 

The second epoch data consist of thirty six 1200 sec images in the F814W filter, and 
provide a baseline of 3 years. The 10 cr limiting magnitude for the point sources in the 
combined second epoch image is AB = 26.7. 

3.1. Data Reduction 

We followed the reduction steps outhned in the HST dither handbook (Koekemoer et 
al. 2002) to reduce the second epoch data. All 36 images were processed with a procedure 
that included the following steps: initial pipeline processing, sky background subtraction. 
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cross-correlation and finding the offsets between images, cosmic-ray rejection, scattered-ligfit 
correction, and final drizzling and combination. 

We obtained the standard pipehne processed images of the HDFS 2001 observations 
from the Space Telescope Science Institute website. After subtracting the sky background 
from each image, cosmic rays were removed so that cross-correlation would not be affected 
by them. Each image is cross correlated with the first image (reference image) in order to 
determine the shifts between them. Most hot pixels, i.e., pixels with elevated dark current, 
can be identified easily with images taken in different dither positions. Nevertheless, some 
hot pixels may escape detection if they fall near an object. We use the static pixel mask 
reference file, used in pipeline calibration, to create a bad pixel mask for our images. 

The cosmic-ray rejection is done in three steps. Each individual input image is registered 
and drizzled to the same output frame. Ten of the 36 images have a higher background and 
suffer from a distinctive cross-shaped pattern due to the shadowing of the HST secondary 
by the support structure of the WFPC2 repeater (Casertano et al. 2000). The 26 images 
without the scattered light problem are combined into a single cosmic-ray-free median image 
for each chip. The median image is then blotted back (reverse drizzling) to the original 
position of each input image. Images with the cross pattern were corrected by applying a 
median filter to the image obtained after subtracting the expected image produced by the 
BLOT task. Inspection of the resulting images showed that the subtraction works well. In 
additon, each image is compared with the blotted images to identify and mask cosmic rays 
and bad pixels. 

The final drizzled images are constructed by applying the shifts and the mask files to 
the individual images, and are combined into a single image for each WFPC2 chip. The 
parameters chosen for the final drizzling were similar to those used for the HDFS 1998 
images. We used a footprint area of 0.6 input pixels (0.5 for the HDFS 1998 images), and 
a pixel scale of 0.4 pixels, resulting in a linear scale of 0.04 " pixel Our experimentation 
with different values of the footprint for the drizzling algorithm showed that image statistics 
are better for a footprint of 0.6 pixels. Because of the relatively small pointing shifts, we 
decided not to combine all images into a single mosaic comprising all four WFPC2 chips. 
The pointing shifts are not large enough to recover part of the sky that is lost between the 
WFPC2 chips. 
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3.2. Object Identification 

Casertano et al. (2000) used the SExtractor package to create source catalogs for the 
HDF South 1998 images. We have created source catalogs for the second epoch data using 
SExtractor, version 2.3, with the same parameters used for the HDF South original images. 
Our main incentive for using SExtractor was its incorporation of weight maps in modulating 
the source detection thresholds. SExtractor follows a standard connected pixel algorithm 
for source identification. We set the source detection threshold to 0.65 and the minimum 
area to 16 drizzled pixels. We convolve the image with a Gaussian (FWHM = 3 pixels) and 
sources with 16 pixels above the detection threshold are included in our catalog. Only those 
objects matching the positions of the objects in the first-epoch data with differences less 
than 1 pixel (0.04") are included in our final catalog. Furthermore, we visually inspected all 
of these sources to avoid any mismatches. We note that the 1 pixel limit is only used for 
the reference compact objects, as we do not impose any limits on the proper motion of the 
point sources. 



4. Proper Motion Measurements 

The original HDF South images were corrected for distortion by Casertano et al. (2000). 
The drizzling procedure that we used to reduce our second epoch data also corrects for 
distortion with standard procedures outlined in the HST dither handbook. Even with these 
distortion corrections, however, some distortion remains (Bedin et al. 2003). Kilic et al. 
(2004) and Pirzkal et al. (2005) showed that a simple quadratic transformation can be used 
to improve the correction for distortion in the WFPC2 and the ACS images, respectively. 

To select reference objects, we use all of the compact objects (isolated, low residuals, 
and not fuzzy) with positional differences less than 1 pixel (0.04") between the two epochs 
to derive a quadratic (third order two dimensional polynomial) transformation for each chip 
(WF2, WF3, and WF4). We rejected deviant points using a 3a rejection algorithm. This 
rejection is required because our reference objects are compact galaxies and centroiding 
errors are larger for galaxies. Even though using a local tranformation for each object can 
increase the accuracy of our proper motion measurements, some of the point sources are 
near the edges of the chip, therefore it is not possible to perform a local transformation for 
each object. Our transformation solutions for individual chips are still better than doing 
a global transformation using a single mosaic image. After mapping the distortions with 
the GEOMAP package, coordinates for the compact objects and the point sources were 
transformed to the second-epoch positions with the GEOXYTRAN routine. 
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Figures 4, 5, and 6 show the observed shifts (scaled by a factor of 50) of each object on 
WF2, WF3, and WF4, respectively. Crosses represent reference compact objects and filled 
circles are point sources identified in the previous section. The transformation error bars 
are shown in the right panel of each figure. It is clear from these figures that the reference 
objects do not show any systematic motion, and several point sources are statistically well 
separated from the reference objects; they are moving with respect to this external reference 
frame. 



4.1. Results 

Our proper motion measurements for the point sources in the HDF South and their 
significances (/u/o") arc given in Table 2. Typical errors in our measurements arc 7-8 mas yr"^. 
Hence, only those objects having proper motions larger than 15 mas yr~^ have significance 
greater than 2(7. Our errors in proper motion measurements are larger than that of Kilic 
et al. (2004) due to the fact that the time difference between the first and second epoch 
images is only 3 years versus 7 years in that study. There are 9 point sources with ^/(t>2 in 
our sample, and four of these stars (HDFS 10081, 10617, 2072, and 895) display significant 
(/x/cr>3) proper motion. 

Figure 7 shows the observed shifts of all point sources along with the Galactic coordi- 
nates. Halo objects rotate around the Galactic center more slowly than disk stars hke the 
Sun. Solar motion relative to the halo objects corresponds to a speed of ~200 km s~^ in 
the direction I = 90°, b = 0° (Mihalas & Binney 1981). Therefore, halo objects are expected 
to lag behind the disk objects in the opposite direction to the rotation direction of the Sun 
(/). A comparison of Figures 4, 5, 6, and 7 shows that even though the reference compact 
objects show a symmetric distribution around the origin, most of the point sources in the 
HDF South are located in the lower left quadrant of the figure; the majority of them lag 
behind the Sun, therefore are hkely to be halo objects. 

Halo stars arc expected to have large proper motions as a result of their high velocities 
relative to the Sun. M&M claimed that the faint blue objects in the HDF South could be 
ancient halo white dwarfs with distances, based on consistency with white dwarf absolute 
magnitudes, less than about 2 kpc from the Sun. Assuming T^an = 200 km s~^ for a halo 
star at a distance < 2 kpc, we would expect to measure a proper motion of > 21 mas yr~^ 
for the faint blue objects. Nine of the faint blue objects identified by M&M (HDFS 1812, 
1827, 1945, 2007, 2178, 441, 1020, 261, 1332) have proper motions in the range 2.26 to 10.27 
mas yr~^, therefore are not likely to be halo white dwarfs. The remaining object, HDFS 
1444, exhibits a significant proper motion, 21.07 ± 7.93 mas yr~^, and is likely to be either 
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a thick disk or a halo object. If the majority of these faint blue objects are not halo white 
dwarfs, what are they? 

5. Spectral Energy Distribution Fitting 

It is always desirable to obtain spectroscopy of star-like objects to distinguish quasars 
from stars, and also to determine the spectral types of stars in order to have a reliable 
estimate of their absolute magnitudes and kinematic properties. Unfortunately, the most 
interesting point sources in deep Hubble images are the faintest, therefore the hardest ones 
for which to obtain spectroscopy. There is only one deep field, the Hubble Ultra Deep Field, 
that is studied spectroscopically to very faint magnitudes (I<27; Pirzkal et al. 2005). Pirzkal 
et al. (2005) found 28 point sources in the Hubble Ultra Deep Field, including 2 quasars 
(and additional two more quasars with I>27). Even though spectroscopy is superior to 
spectral energy distribution fitting using photometry, Kilic, von Hippel, & Winget (2005) 
could identify the quasars and the stars in the HUDF from their photometry. Using the 
IRAF package CALCPHOT, which is designed for simulating the HST observations, Kilic, 
von Hippel, Sz Winget (2005) showed that spectral energy distribution (SED) fitting can be 
used to identify quasars, and to assign spectral types to stars. Their classifications for the 
stars agreed reasonably well with the Pirzkal et al. (2005) spectroscopic classifications, and 
they were also able to predict the redshifts of two of the quasars correctly, whereas the other 
two were predicted to be at lower redshifts than the spectroscopic redshift measurements. 

Here we adopt the same SED fitting procedures used by Kilic, von Hippel, & Winget 
(2005). We use Pickles (1998) stellar templates and the composite-quasar spectrum from 
the Sloan Digital Sky Survey (Vanden Berk et al. 2001) to simulate SEDs in U300, B450, 
Veoe, and I814 using the CALCPHOT task. Our simulations include stellar templates from 
05 to M6 dwarf stars, and the composite quasar spectrum is used to simulate the colors for 
quasars up to 2; = 4.8. Note that we do not simulate the colors of the broad range of galaxy 
types, as almost every galaxy with sufficient signal-to-noise is resolved by HST (e.g. 2 kpc 
at z = 0.5 projects to 0.33" for Qo = 0.3, Ao = 0.7, and Ho = 70). 

The observed magnitudes for the faint blue objects are converted to F^, normalized at 
Veoe, and compared with our simulated SEDs. We assign spectral types to each object using 
a minimization technique. Each photometry point is weighted according to its errorbar, 
but we also tried giving equal weights to each photometric band to explore possible fits. 
Figure 8 shows our best-fit solutions to the point sources in the HDF South. Observed 
fluxes are shown as filled circles, the solid lines show the best fit stellar or quasar templates 
when each photometric band is weighted according to its errorbar, and the dashed lines 
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show the best fit templates assuming equal weights for all bands. The object number and 
the assigned spectral types are given in the top left corner of each figure. Error bars for 
the observed fluxes are so small that they can only be seen for the fainter objects. We note 
that the U and B photometry for HDFS 2615 is probably wrong, as it is barely detected 
in the U and B images, yet its photometry shows excesses in the U and the B filters, also 
apparent in Figure 13. Even though our SED fitting algorithm gives a range of spectral 
types (09-M2) for this object, it is more hkely to be an M type star based on more accurate 
V and / photometry. 

Figure 8 shows that there are 4 point-like sources with SEDs better explained by quasars 
than stars. HDFS 10151 {z = 2.0), 1945 {z = 0.2), 2007 {z = 4.7), and 2178 {z = 4.0) have 
colors more consistent with being quasars. Interestingly enough, three of these objects were 
classified as possible halo white dwarf candidates by M&M. These four objects have proper 
motions in the range 2.56 - 4.91 ± 8 mas yr~^ (see Table 2), and are not moving (within 
the errors); they demonstrate that our proper motion measurements are reliable. 

6. Discussion 

We identified 4 quasars and 42 stars, including 36 KO or later type stars, with our SED 
fitting technique. We adopt the absolute magnitude (My) for each spectral type from Pickles 
(1998; Table 2). We use the My for each object to simulate the absolute magnitudes in the 
WFPC2 F606W filter (Mpeoew)- We calculated photometric distances for all of the objects 
in our sample using the simulated Mf606W and the observed Veoe magnitudes. Combining our 
proper motion measurements with photometric distances, we are able to derive tangential 
velocities for the stars. The photometric errors for the point sources in the HDF South are 
small, hence they cause only relatively small errors in the estimated distances. On the other 
hand, if the assigned spectral types are wrong by 1 index (for example M3 instead of M2), 
then the absolute magnitudes could be wrong by as much as 1 magnitude. We also note 
that our stellar templates have approximately solar metallicity. Therefore, the distances to 
the metal poor halo objects, which will be intrinsically fainter for the same spectral type, 
may be over-estimated by our SED fitting method. For example, a metal poor GO dwarf 
{[Fe/H] = —0.8) would be 0.5 magnitude fainter than a solar metallicity GO star (Pickles 
1998). Hence, the distance to a metal poor GO type star would be over-estimated by 26%. 

In order to determine the effect of different metallicites, we use synthetic spectra from 
a PHOENIX model atmosphere grid (Brott & Hauschildt 2005) for stars with 2000 K < 
Tcs < 10000 K, [Z/H]= 0.0, -0.5,-2.0, and log(g)=4.5 to simulate photometric colors in 
the HST bands. We found that using different metallicities changes the best fit T^s by 
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50 - 450 K for GO and later type stars, and the spectral types obtained from metal poor 
atmosphere models are usually 1-2 later in spectral types, e.g. M2 instead of MO, than 
the ones obtained from the models with solar metallicity. Wc also used Kurucz model 
atmospheres (Kurucz 1995) with 3500 K < T^s < 6000 K to test the metallicity effect, and 
found similar results. Using Girardi et al. (2002) theoretical isochrones, we estimate My and 
Mf606W for different metallicities for the point sources in our sample. Spectral types from the 
SED fitting procedure (using the Pickles library), estimated absolute magnitudes, derived 
distances, and tangential velocities are given in Table 3. The ranges of absolute magnitudes, 
distances, and tangential velocities show the effect of using models with different metallicities. 
Despite these potential systematic errors, the total information we have available for the point 
sources (the morphological information, proper motion measurements, and the SED fitting 
results) is sufficient to determine their nature. 



6.1. Late Type Stars in the Galaxy 

Figure 9 shows the histogram of the number of stars observed at a given distance in 
the HDF South along with the predictions from the star count models (Reid & Majewski 
1993; dashed line) and the observed distribution of stars in the HUDF (Pirzkal et al. 2005; 
solid line). The top panel shows the distribution of stars if we use Pickles main sequence 
library, and the bottom panel shows the same distribution if we use synthetic spectra with 
[Z/H]=-2.0 (halo-like metallicity). Both panels show that there is a deficit of stars in the 
5~15 kpc range in the HDF South compared to the star count models. 

If we assume that all of the stars in the HDF South have solar metallicity, then we 
see an excess number of stars at ~35 kpc, and 4 more stars (HDFS 441, 1020, 2596, and 
261) with estimated distances larger than 50 kpc, which is not expected from the star count 
models, nor would the distances be consistent with the metallicities for Galactic stars. Most 
of the objects observed at ~35 kpc are K type stars, whereas the four objects with distances 
> 50 kpc fit late type star SEDs fairly well, and they have essentially the same spectral type 
(K7-M0). If we assume that their classification as point sources is reliable, our SED fitting 
procedure works well to identify quasars and stars, and that they have solar metallicity, then 
we could claim that we discovered two new populations of stars; a cluster of stars at 35 kpc 
and several other stars at > 50 kpc. 

HDF South and HUDF have similar Galactic latitudes (-49.21 and -54.39, respectively), 
therefore they should have similar distributions of Galactic objects. We do not find any stars 
with distances larger than 50 kpc in the HDF North or the HUDF. Even though HDF South 
(/ =328.25, 6 =-49.21) is about 25° away from the center of the Small Magellanic Cloud 
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(SMC; I =302.80, b =-44.31), the Magellanic Stream has over dense regions near the HDF 
South (Ix 10^^ atoms cm~^; see the HI surface density maps of Mathcwson & Ford 1984). 
The average distance to the 4 objects with d > 50 kpc is 66.8 ± 13.4 kpc (assuming solar 
metallicity), whereas the distance to the SMC is measured to be 60.6 ± 3.8 kpc (Hilditch, 
Howarth, & Harries 2005). Can these objects be members of the SMC? Assuming that they 
have metallicities similar to the SMC {[Fe/H] = —0.7; Lennon 1999), we estimate their 
average distance to be 17.8 ± 3.4 kpc, therefore they are not consistent with being in the 
SMC. 

We expect the majority of the stars in the HDF South to be halo stars (from the 
star count models and the velocity distribution of Figure 7). Therefore, the bottom panel 
of Figure 9 is hkely to represent the real distribution of stars better than the top panel. A 
comparison of the distances, tangential velocities, and photometric colors show that all of the 
42 stars except three (HDFS 1444, 895, and 2488) are consistent with being dwarf /subdwarf 
stars in the Galaxy. Seven of these 42 stars (HDFS 191, 1576, 1922, 2072, 10081, 10326, and 
10617) show significant proper motion {/i/a > 2) and have spectral types M2-M5, therefore 
are probable halo M dwarfs. Large errors in our proper motion measurements prevent us 
from classifying the kinematic properties of the other stars, nevertheless, they are most likely 
to be GO and later type dwarfs in the thick disk or halo of the Galaxy. The three stars (HDFS 
1444, 895, and 2488) with estimated distances larger than 90 kpc (both for the metal rich 
and the metal poor case) are discussed in the next section. 



6.2. White Dwarfs 

HDFS 1444, 895, and 2488 would have to be at very large distances (> 90 kpc) if 
they were main-sequence stars of any metallicity. On the other hand, as white dwarfs, 
they would be at more reasonable distances. In order to find the temperatures of these 
objects, we simulated the colors for blackbody SEDs with temperatures in the range 3000 - 
80000 K. Figure 10 shows the best fitting blackbody SEDs (solid fines) for HDFS 1444 (top 
panel), 895 (middle panel), and HDFS 2488 (bottom panel). We calculate the blackbody 
temperatures for these objects to be 10547 K, 6096 K, and 10463 K, respectively. Wc have 
also used DA white dwarf models (D. Saumon and D. Kocster, private communication) 
to simulate colors for white dwarfs with log g = 8 and 3000 K < T^s < 20000 K. Cool 
white dwarfs show depressed infrared fluxes due to the effects of collision induced absorption 
(CIA) due to molecular hydrogen (Hansen 1998, Saumon & Jacobson 1999). The pure H 
white dwarf models that we used include the CIA opacities, therefore, we are able to compare 
the spectral energy distributions of young and old white dwarfs simultaneously and find the 
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best-fit solution for our white dwarf candidates. Assuming that HDFS 1444, 895, and 2488 
are pure-H atmosphere white dwarfs, we estimate the temperatures of these objects to be 
10681 K, 5882 K, and 11000 K, respectively. Our best fitting DA white dwarf SEDs are 
shown as dotted lines in Figure 10. 

Using our best fit DA white dwarf atmosphere solutions, we estimate the absolute 
magnitudes for our white dwarf candidates using the tables from Bergeron et al. (1995). We 
use the white dwarf models to predict Mf606W! and therefore to calculate the distances and 
tangential velocities for these three objects. Estimated distances and kinematic properties 
of our white dwarf candidates are given in Tabic 3. HDFS 1444 displays a significant proper 
motion, 21.07 ± 7.93 mas yr and is consistent with being a thick disk or halo object at 
1.5 kpc with Vts,n = 148 ± 56 km s~^. Likewise, HDFS 895 displays a proper motion of 34.9 
± 8.0 mas yr and is more likely to be a halo white dwarf at 2.1 kpc with T^an = 346 ± 
79 km s~^. HDFS 2488 does not display any significant proper motion, and its classification 
as a point source is questionable (see §2.1), nevertheless, if it is a star, then it would have 
to be a halo white dwarf at 9 kpc. 

Figure 11 shows the U vs. V velocity diagram for the stars (triangles; assuming [M/H]=- 
2.0) and the likely white dwarfs (circles) in the HDF South. We use the results of Chiba 
& Beers (2000) for the expected velocity distribution of halo (solid line) and thick disk 
(dashed line) objects. It is apparent from this figure that we have several thick disk stars 
in our sample. Nevertheless, the majority of the stars are more likely to be in the Galactic 
halo. One of the white dwarf candidates, HDFS 1444, may well be a thick disk white dwarf, 
whereas the other two candidates, HDFS 895 and 2488, are more likely to be halo white 
dwarfs. 



6.3. Mendez &; Minniti Faint Blue Objects 

M&M identified 10 faint blue objects in the HDF South. We classified two of these 10 
objects (HDFS 1812 and 1827) as resolved (see Figure 2). We measure proper motions of 
2.26 ± 7.93 and 3.25 ± 6.98 mas yr^^ for these two objects, respectively. One of the faint 
blue objects (HDFS 1332) is fainter than Isu — 27 (therefore not included in our analysis), 
and has a proper motion of 9.31 ± 7.93 mas yr~^. Six of the faint blue objects (HDFS 1945, 
2007, 2178, 441, 1020, 261) have proper motions in the range 2.56 to 10.27 mas yr'^ Our 
SED fitting analysis showed that HDFS 1945, 2007, and 2178 have colors more consistent 
with being quasars than stars. In addition, we classify HDFS 441, 1020, and 261 as metal 
poor stars in the halo of the Galaxy. Therefore, only one of the faint blue objects identified by 
M&M, HDFS 1444, plus two more white dwarf candidates identified in our analysis (HDFS 
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895 and 2488) are consistent with being white dwarfs. 

We use Reid & Majewski (1993) star count models and our own calculations based on 
Gilmore et al. (1989; see Kihc et al. 2004 for a complete discussion) to predict the number 
of stars and white dwarfs expected in the HDF South. We expect to find 45-52 stars and 
0.66-2.31 white dwarfs, including 0.24-0.50 disk white dwarfs in the HDF South. The star 
count models mildly over-predict the observed number of stars. The observed population of 
2-3 white dwarfs is consistent with the standard Galactic models. 



7. Conclusion 

Superb resolution of the HST WFPC2 camera helped us to ehminate almost all of the 
large number of galaxies present in the HDF South. Using accurate SED fitting procedures 

(to eliminate quasars) and obtaining proper motion measurements with a 3 year baseline 
(to identify fast moving-halo objects) enabled us to classify the point sources in the HDF 
South. We identified 4 quasars (consistent with zero proper motion) and 42 stars. Three of 
these stars (HDFS 1444, 895, and 2488), if on the main sequence, are too distant to be in the 
Galaxy, and are best explained as white dwarf stars. Their kinematic properties show that 
HDFS 1444 is probably a thick disk object, whereas the other two are more fikely to be in 
the halo of the Galaxy. Third epoch data on the HDF South would be useful to place better 
constraints on the proper motions and the kinematic memberships of the point sources. 

Figure 12 and 13 show the location of the stars, the quasars, and the white dwarf 
candidates in color-magnitude and color-color diagrams. A comparison of Figure 12 with 
Figure 3 of M&M shows that a detailed analysis of the faint blue objects is needed to 
classify their nature-important additional constraints can be derived from SED fitting and 
from proper motions, even with a baseline of only 3 years. The number of blue extra-galactic 
sources is enormous at these magnitudes compared to the number of stellar objects. Our 
analysis shows that only two of the faint blue objects (HDFS 1444 and 895) show significant 
proper motion and are consistent with being disk/halo white dwarfs. None of the other faint 
blue objects exhibit significant proper motion, and therefore they are highly unlikely to be 
halo white dwarfs. Even if we misclassified a few blue extra-galactic objects as unresolved, 
this would only decrease the observed number of stars and white dwarfs in the HDF South, 
and therefore, strengthen our conclusion that the majority of the faint blue objects are not 
halo white dwarfs. The observed population of ~2 white dwarf candidates (/ < 27) in the 
HDF North (Kilic et al. 2004), HDF South (this study), and HUDF (KiUc, von Hippel, & 
Winget 2005) imply that the faint blue objects, and especially the non-observed population 
of halo white dwarfs, are highly unlikely to solve the dark matter problem. 
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Table 1. Point Sources in the Hubble Deep Field South 
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Tablc 1 — Continued 
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Table 2. Proper Motions 
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316 


wf4 


-0.10 


-0.38 


5.29 


8.00 


0.66 


1426 


wf2 


-0.91 


-0.45 


13.50 


7.93 


1.70 


652 


wf4 


-0.31 


-0.36 


6.32 


8.00 


0.79 


191 


wf4 


-0.98 


-0.78 


16.73 


8.00 


2.09 


1302 


wf3 


-0.34 


-0.20 


5.22 


7.04 


0.74 


895 


wf4 


-0.34 


-2.60 


34.90 


8.00 


4.36 


2615 


wf2 


0.32 


-0.48 


7.62 


7.93 


0.96 


1945 


wf3 


0.30 


0.21 


4.91 


7.04 


0.70 


2007 


wf3 


-0.15 


-0.17 


3.03 


7.04 


0.43 


2596 


wf3 


-0.03 


-0.14 


1.90 


7.04 


0.27 


441 


wf4 


-0.75 


-0.19 


10.27 


8.00 


1.28 


1020 


wf4 


-0.44 


-0.57 


9.63 


8.00 


1.20 


2178 


wf2 


0.14 


-0.13 


2.56 


7.93 


0.32 


261 


wf4 


0.35 


-0.27 


5.94 


8.00 


0.74 


1306 


wf2 


-0.54 


0.09 


7.34 


7.93 


0.93 


2488 


wf3 


-0.38 


-0.17 


5.60 


7.04 


0.80 



Note. — HDFS 15, 1010, 1187, 1209, 2257, and 10017 are out of the field of 
view in the second epoch data. 
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Table 3. Spectral Types, Photometric Distances, and Tangential Velocities 



Object 


Type 




d(kpc) 


14„„(kni/s) 


Wood 


1 TV flT\ 
Kl-MO 


C O OCT 

5.73-8.25 


1 C AO 

1.5-4.8 


O A O^JA 

84-269 


10692 


A3-4 


1.42-4.97 


7.1-36.4 


1 OA f\ A Ci 

189-948 


10017 


K5-6 


'TOO 1 r\ p" o 

7.33-10.58 


0.9-3.9 




iUooo 


A /TO 

M2 


HIT 11 nn 

y.iY-ii.uy 


c\ o 1 n 

u.o-i.y 


36-87 


iUioi 










10081 


IV JVC 

JV15 


10 10 10/1/^ 


0.9-1.0 


1 1 O A O A 

118-134 


1576 


M3 


10.19-12.56 


0.7-2.2 


75-222 


2257 


M3 


10.19-12.56 


1.2-3.4 




iUoiv 


Mo 


1010 1 ^ AC 

12.18—12.46 


1.7-1.9 


1 OO 01 A 

188—214 


8d 


G3 


4.65-6.44 


12.0-26.4 


o /) r\ COT 

240—527 


15 


M3 


in in locr/^ 

10.19 12.56 


1.4 4.2 




2041 


M5 


12.18-12.46 


2.2-2.5 


67-76 


10326 


TV /TO 


9.17-11.09 


2.3-5.5 


A 1 P* AO 

216-522 


7Ui 


A /TO C 

Mz.o 


n 11 Tn 
y.oo-ii.7U 


on CO 


1/1/1 o/:! c 

144-365 


1 noo 
1922 


A /TO 


n 1 T 11 on 

y. 17-11. uy 


Z.O-D.U 


206-498 


1209 


Mb 


10.86-13.56 


1.7-6.0 




1257 


JV15 


1010 10 AO 

12.18-12.46 


2.6-3.0 


119-136 


431 


Kl 


5.73-7.24 


15.1-30.3 


AAA A C\f\ 

200-400 


1187 


Ml 


O fCX 11 1/1 

8.59-11.14 


3.4-10.5 




2469 


K.1 


5.73-7.21 


1 C A 1 C 

15.9-31.5 


A Af\ OT /I 

442-874 


323 


Kl-3 


p- T /to 

5.73-7.43 


15.7-34.4 


A t~l A A /t A 

434-949 


1 AO 

108 


TV /ri 

Ml 


D /^n 11 1 /t 

8.09-11.14 


O T 11/^ 

3.7-11.6 


OO T TO 1 

237-731 


1444 


B8-9 


-1.39 - -0.09 


379.2-690.0 


OTOTr *^OAA1 

37875-68921 


^ A A A 
14:4:4: 


WU 


11.95 


1.482 


1 /I O _L Cf? 

148 ± 56 


2364 


M5 


12.18-12.46 


3.7-4.2 


48-54 


1724 


K3 


6.58-7.68 


15.2-25.3 


324-537 


1386 


K6 


7.59-10.64 


4.4-18.0 


231-942 


933 


M3 


10.19-12.00 


3.6-8.2 


109-251 


2072 


M5 


12.18-12.46 


4.0-4.5 


399-454 


1331 


M2 


9.17-11.09 


5.2-12.6 


194-469 


316 


K7 


7.86-11.15 


4.8-22.0 


121-552 


1010 


K3 


6.58 7.41 


23.1-33.9 




1426 


MO 


8.25 11.13 


6.1-22.8 


388 1460 


652 


K6 


7.59-10.50 


7.9-30.1 


236-901 


191 


M5 


12.18-12.46 


7.4-8.5 


589-670 


1302 


MO 


8.25-11.13 


7.8-29.3 


193-725 


895 


GO 


4.18-5.98 


93.5-208.4 


15469-34471 
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Table 3 — Continued 



Object 


Type 




d(kpc) 


14a„(km/s) 


895 


WD 


14.17 


2.094 


346 ± 79 


2615 


09-M2 


-4.32-11.09 


14.6-17600 


527-636472 


1945 


QSO 


z=0.2 






2007 


QSO 


z=4.7 






2596 


MO 


8.25-11.13 


14.2-53.3 


127-479 


441 


MO 


8.25 11.13 


15.7 59.3 


766-2886 


1020 


K7 


7.86-11.15 


15.6-71.1 


713-3244 


2178 


QSO 


z=4.0 






261 


MO 


8.25-11.13 


22.2-83.5 


623-2349 


1306 


M3 


10.19-12.18 


18.1 45.4 


631-1577 


2488 


AO-1 


0.71-1.00 


1350-1540 


35803-40918 


2488 


WD 


11.86 


9.070 


241 ± 303 
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Fig. 1. — The distribution of R50 (half hght radius) for the HDF South objects. The sohd 
hne separates the unresolved objects from the resolved objects. 
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Fig. 2. — Point spread function distributions of the objects with ^ + / > 21, R50 < 3 pixels 
and stellarity index > 0.85. The object number, R50, FWHM, and stellarity index of each 
object are shown in the top right corner of each paneL The HDFS 10617 PSF is shown in 
each panel for comparison. 
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Fig. 2. — cont. 
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Fig. 2. — cont. 
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Fig. 3. — The relation between R50 and FWHM for the HDFS objects in different magnitude 
ranges. Objects with stellarity index >0.85 arc shown as filled triangles, whereas the open 
triangles show the rest of the objects in the HDFS catalog. Filled circles represent unresolved 
objects identified from their empirical PSF distributions. Point sources can be identified 
easily down to V + / = 27, and the star/galaxy sequences begin to blur for ^ + / > 27. 
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Fig. 4. — Left panel: observed shifts (scaled by a factor of 50) of point sources (filled circles) 
and surrounding reference compact objects in the HDFS WF2 chip between Epoch 1 and 
Epoch 2 images; Right panel: Differences between the first-epoch coordinates (Xi, Yi) and 
the transformed second-epoch coordinates (X2, Y2) for the same objects. The expected 
scatter caused by the transformation itself is represented by the error bar in the lower right 
corner of the figure. The axes units are HDFS pixels (0.04"/pixel). 
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Fig. 5.- 



Same as Figure 1, but for the objects detected in WF3 
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Fig. 6. — Same as Figure 1, but for the objects detected in WF4 
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Fig. 7. — Observed shifts of all point sources in the HDF South. The directions of increasing 
Galactic longitude (/) and latitude (6) and the South Galactic Pole (SGP) are also shown. 
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Fig. 8. — Spectral energy distribution and best fitting template for the point sources in the 
HDF South. 
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Fig. 8.— 



cont. 
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Fig. 8.— 



cont. 
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Fig. 9. — Histogram of the number of stars observed at a given distance in the HDF South 
(shaded histogram) assuming solar (top panel) or metal poor composition (bottom panel). 
The histogram of the number of stars observed in the HUDF (solid line) and the predictions 
from the star count models (dashed-dotted line) are also shown. 
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Fig. 10. — Fits to the spectral energy distributions of HDFS 1444, 895, and 2488. Best 
fitting blackbody models are shown as solid lines, whereas bestfitting pure-H atmosphere 
white dwarf models are shown as dotted lines. The first temperature estimate in the label 
is for the blackbody and the latter one is for the pure-H white dwarf models. 
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Fig. 11. — U vs. V diagram for the stars (filled triangles) and the hkely white dwarfs (filled 
circles) in the HDF South. The 2(7 velocity ellipse of the thick disk objects is shown (dashed 
line) along with the la ellipse of the halo population (solid line). The errors in the tangential 
velocities are larger for distant objects, therefore these objects have also large errorbars in 
the U vs. V diagram. 
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Fig. 12. — The location of the 46 point sources in a color magnitude diagram. Filled circles 
and triangles represent stars and quasars identified from spectral energy distribution fitting, 
respectively. White dwarf candidates are marked with open circles, whereas the smaller dots 
show the rest of the objects in the HDF South catalog. 
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Fig. 13. — Vqoq — /775 vs. ^435 — VeoQ color-color diagram for the stars (circles), the quasars 
(triangles) and the white dwarf candidates (open circles) in the HDF South. Predicted 
tracks for main sequence stars (05-M6; dotted line), pure-H white dwarfs {log g=8, T^s — 
60000-3000K; dashed line), and quasars (sohd line) are also shown. 



